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We consider the transport of rigid objects with internal structure in a flashing ratchet potential 
by investigating the overdamped behavior of a rod-like chain of evenly spaced point particles. In 
ID, analytical arguments show that the velocity can reverse direction multiple times in response to 
changing the size of the chain or the temperature of the heat bath. The physical reason is that the 
effective potential experienced by the mechanically coupled objects can have a different symmetry 
than that of individual objects. All analytical predictions are confirmed by Brownian dynamics 
simulations. These results may provide a route to simple, coarse-grained models of molecular motor 
transport that incorporate an object's size and rotational degrees of freedom into the mechanism of 
transport. 



I. INTRODUCTION 

A number of theoretical and experimental studies in 
recent years have addressed directed transport of diffu- 
sive particles in spatially periodic systems in the absence 
of net forces, which occurs when spatial or temporal in- 
version symmetry is broken while the system is kept away 
from thermal equilibrium P, |3| . Such systems are called 
ratchets or Brownian motors [2| , and their study has both 
fundamental and practical motivations, including appli- 
cations in biology and nanotechnology jj] . There are var- 
ious mechanisms by which the system can be maintained 
out of thermal equilibrium, including a time-dependent 
force (often called a "rocking" ratchet) 0, 0, Q and a 
time-dependent potential barrier (a "flashing" ratchet) 

m 

While many ratchet studies have dealt with the asym- 
metric pumping of individual, pointlike particles 0, 0, 
ITol ITU , motors with internal structure have recently at- 
tracted interest. One reason is that such studies may 
provide a route to models of the linear transport of bi- 
ological molecular motors. A number of studies have 
demonstrated behavior of mechanically coupled particles 
which is qualitati vely different than that of a s ingle pa rti- 
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|2^ l28l l2^ l30| ■ for example: In contrast to the behavior 
of an individual particle, two harmonically coupled par- 
ticles in a flashing ratchet undergo directed motion in 
the absence of thermal fluctuations 0, ^3 ■ Fo^' a flash- 
ing ratchet, two harmonically coupled particles have a 
smaller velocity than a single particle, while for a rock- 
ing ratchet, the coupled particles have a greater velocity 
|l4| . Two rigidly coupled particles in the presence of 
a stochastically rocked ratchet potential undergo trans- 
port that reverses direction as a function of the dimer 
size ^1 Large collections of mechanically coupled 
particles in a flashing ratchet can undergo spontaneous 
oscillations, which may be relevant to the collective mo- 
tion of molecular motors in muscles ^3 ■ 

In this study, we wish to establish a conceptual un- 



derstanding of the role of coupled motion in a flashing 
ratchet by considering the most simple possible form of 
mechanical coupling between particles: a chain of rigidly 
connected point particles. The chains are exposed to 
Gaussian white noise in the overdamped limit, and driven 
by a periodically modulated, asymmetric, spatially peri- 
odic potential (Fig.^. The periodically flashing ratchet 
scheme is chosen because it is simple enough to allow an 
intuitive understanding and analytical prediction of the 
behavior of coupled particles. 

Using analytical arguments, we show that, for a ID 
system, the average velocity can reverse direction in re- 
sponse to changing the size of the chain or the temper- 
ature of the heat bath without changing the symmetry 
of the applied potential. However, when the chain is al- 
lowed to rotate freely in 3D, or when its length is much 
less than the spatial period of the ratchet potential, ve- 
locity reversal is no longer observed, and the qualitative 
behavior of single particle motion is recovered. The be- 
havior observed for mechanically coupled particles in a 
flashing ratchet can be understood in terms of the chain's 
center-of-mass effective potential, which can have differ- 
ent symmetry than the potential felt by an individual 
particle. All analytical predictions are confirmed by nu- 
merical simulations. 

The layout of this paper is as follows: In sectional we 
introduce the Brownian dynamics (BD) model used in 
our simulations. In section UllI we discuss the motion of 
chains of two or more evenly spaced point particles con- 
strained to ID motion in the limit of perfect confinement 
{kT < Knax)- In sectionllVl we discuss how the direction 
of ratchet velocity for two coupled particles moving in 
ID depends on the temperature of the surrounding heat 
bath. The behavior of rigidly connected particles in 3D is 
discussed in the section Ivl In the conclusion, we discuss 
a possible experimental realization of this system, which 
may be used to directionally separate polymer segments 
of different length. We also discuss the possible biologi- 
cal relevance of our results, namely the observation that 
complex molecular motors with only small differences in 
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structure can move in opposite directions. 
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FIG. 1: The applied ratchet potential V{x), shown above, is 
characterized by periodic length, L, height, Vmax and asym- 
metry, a. 



II. THE MODEL 

We consider mechanically coupled particles periodi- 
cally subjected to a piecewise linear 'sawtooth' potential, 
V{x) (Fig.^, characterized by periodic length, L, height, 
VInax, and asymmetry, a. The particles are alternatingly 
subjected to V{x) for a time period ton, and allowed to 
diffuse freely for a time period ioff, thus asymmetrically 
harnessing Brownian motion to produce net transport 
which, for a < 1/2, is in the +x direction in Fig. H for 
non-interacting particles (see section Hill and 2, 8, 12]). 

The simulations in this study are based on the follow- 
ing scheme: The coupled particles are modeled as a chain 
of point-like beads, separated by a distance d, defined by 
a repulsive Lennard-Jones interaction 

[ : rij > 26(7 

(1) 

where is the separation between beads, and a finite ex- 
tensible nonlinear elastic (FENE) potential between ad- 
jacent beads: 



U{n,) = --kpRlln l^l - . (2) 
The equations of motion of individual beads are given 

by 

mr\ = -7f,ri + ^{t) - VViNrin) - WsxTit, n), (3) 

where ^{t) is a randomly fluctuating Gaussian white 
noise term with zero mean and correlation (C(t)C(s)) — 
27f,kT(5(t-s), 76 is the drag coefficient of a bead, k is the 
Boltzmann constant and T is the temperature of the heat 
bath. The term ViNri^i) represents the potential of the 



bead due to its interaction with the other beads in the 
chain, and VExxitj^i) is the external ratchet potential. 
The equations of motion are integrated using a Brownian 
dynamics algorithm. The mass, m, of each bead is set to 

unity and we use a, e and r = \J~^^ as scaled units of 
length, energy and time. 

Unless otherwise specified, L — 5cr, and i?o is cho- 
sen such that the distance between adjacent beads is 
d = 0.97f7 for kT = e. In sections ImH and UTTbI the 
FENE parameter i?o is altered to adjust the value of d. 
In section IIVI the temperature T is varied, but i?o is 
adjusted accordingly, so that d is the same for different 
values of the temperature. We set ion = ^off = 20, 
because this gives a chain of several beads enough time 
to localize during ion and to diffuse to an adjacent well 
during toff- 

In order to isolate the role of an object's geometry in 
ratchet transport from the effect of varying the diffusion 
constant, we give all chains the same total drag coeffi- 
cient. The drag coefficient for one bead is chosen to be 
7h — 1/N, where N is the number of beads in the chain. 
Since hydrodynamic effects have not been included, the 
chain's total drag coefficient is 7t = N^i, = 1 31]. 

When the particles are constrained to motion along the 
X-axis, this model describes a rigid rod made up of evenly 
spaced particles. When individual particles are allowed 
to move in three dimensions, this system corresponds to 
the Rouse model of a polymer, with equilibrium config- 
urations described by a self-avoiding walk [TsL l3^ . This 
modeling scheme is used here for rigid rods in order to 
provide continuity with our investigation of the role of 
flexibility in flashing ratchet transport [l8| . 

III. EFFECTS OF COUPLING IN LOW 
TEMPERATURE REGIME 

Before we describe how mechanically coupling the mo- 
tion of point particles affects their ratchet velocity, we 
briefly review the motion of a single particle (a monomer) 
in a flashing ratchet 0, 0, ^3 ■ ^^e low-temperature 
limit {kT <C Knax), a monomer will have positive veloc- 
ity when a < 1/2 for the following reason: During ton, 
the monomer localizes at the minimum of the sawtooth 
potential: Xmin = (1 — ct)I/, where a: = corresponds 
to the beginning of the potential well. The minimum 
distance the monomer must diffuse during toS in order 
to localize in the adjacent well in the -l-a; direction dur- 
ing ton is Aa;4. = aL. The minimum diffusion distance 
to localize in the adjacent well in the —x direction is 
Ax- = (1 — a)L. The ratchet velocity will be positive 
if Aa;+ < Ax-. As the asymmetry, a, is increased from 
zero, there is a reversal from positive to negative velocity 
at a critical asymmetry, ac = 1/2, given by the condition: 
Aa:+ = Ax-. 

For a ID chain of particles in the limit kT <C Knax, the 
characteristic diffusion distances, Ax_|_ and Ax^, can be 
calculated for the chain's center of mass (CM) by consid- 
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ering the shape of the effective potential, U{xcm), given 
by 



1 ^ 

U{xcm) = j;^J2^ix.O. 



(4) 



The center of mass of a chain of coupled particles will 
localize exactly at a local minimum, Xmin, of its effec- 
tive potential. The distance Ax+ (Ax_) is given by the 
distance between x^i^ and the closest absolute effective 
potential maximum in the +x {—x) direction (Fig. 
The net ratchet velocity is positive (negative) for a < Uc 
(a > Qfc), respectively, where Ofc is again given by the 
condition Aa;+ = Aa;_. For simplicity, we limit our dis- 
cussion to chains with total length less than a spatial pe- 
riod of the ratchet {{N — l)d < L). The velocity of any 
object is antisymmetric about a — 1/2, so only values of 
a in the interval [0; 0.5] are discussed in this study. 



A. Dimer 

We begin to explore the role of coupled particle motion 
in a flashing ratchet by considering two coupled particles 
(a dimer) separated by a distance d, constrained to ID 
motion. The system behaves differently depending on 
whether d < aL or d > aL, and we will briefly discuss 
each of these cases. 

When d < aL, the dimer localizes to the center-of- 
mass position a;inin = (1 — a)L — d/2 during ton because, 
at this position, one bead experiences the shallow slope of 
V{x) and the other sits at the minimum of the potential 
well. Any displacement in the +x or —x direction would 
lead to a restoring total force, and so Xmin corresponds 
to a minimum of U{xcm) (Fig.EJa)). 

For a dimer to localize in the adjacent well in the +x 
direction after a ioff period, it must diffuse until xqm > 
L + d/2, such that both beads are in the next well and 
the dimer experiences a total force Fqm > when the 
potential turns on. Therefore, Ax^ = L + d/2 — Xmin ~ 
aL + d. 

Likewise, for a dimer to localize in the adjacent well 
in the —x direction after a ioff period, it must diffuse 
until xcM < d/2, such that one bead is exposed to the 
steep side of the adjacent well in the —x direction, and 
^CM < at the beginning of ton- Consequently, Ax- = 
Xmin — d/2 = (1 — a)L — d, and the condition, Axj^ = 
Ax^, yields the critical asymmetry 
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2 L 



(5) 



Since Eq. © was obtained by assuming d < aL, it 
holds for d/L <a^^\/2- d/L, or d/L < 1/4. 

When d > aL, the effective potential U{xqm) has two 
local minima: Xmini = (1 ~ C()L — d/2 and Xmin2 = 
(1 — a)L + d/2, where Xmmi corresponds to localization 



inside a V{x) potential well, and Xmin2 corresponds to a 
dimer straddling a V{x) potential maximum (Fig. I^Jb)). 
The average position during ton depends on the relative 
probability for the object to be trapped in each of these 
minima, which is a function of a, d/L and toff- If the 
dimer has a probability $ (a, j-, toff) of being trapped in 
the position Xmini at the end of tos , then the average lo- 

Cahzation position is: (Xmin) = ^Xminl + (1 - ^)Xmin2 = 

(l-a)L + (l-24>)d/2. 

In this case, the characteristic diffusion distances are 
Ax-|_ = L + d/2 — (xmin) and Aa;_ = (xmin) — d/2, and 
the condition, Aa;+ = Ax^, yields 



^ ^ f d \ d 
ttc = 2 - * I "c, ioff 1 (6) 



Note that for d < aL, where $ = 1, Eq. © simplifies to 
Eq. ©. 

These analytical predictions are confirmed by BD sim- 
ulations. Figure Efa) shows simulation results for the 
time-averaged velocity of dimers with different d, for 
kT/Vraax = 1/50, demonstrating that the ratchet ve- 
locity reverses as a function of a. Figure El^b) confirms 
the prediction of Eq. that ac decreases linearly with 
d/L, for d/L < 1/4, and is no longer linear with d/L for 
d/L > 1/4 (Eq. ©). 



B. Three or more coupled particles 

We now extend our discussion to chains of three or 
more particles moving in ID. In particular, we make an- 
alytical arguments for when a velocity reversal with a is 
expected for three coupled particles (a trimer), and sup- 
port these predictions with BD simulations. Analytical 
predictions become increasingly complicated for iV > 3, 
and we simulate chains with A'^ = 4 and A^ = 5 to make 
general observations about ID coupled motion in a flash- 
ing ratchet. 

For A^ coupled particles, each separated by a distance 
d, the number and locations of minima in the effective 
potential depend on the following factors: (1) Is the po- 
tential asymmetric enough that one particle on the short 
(steep) side will experience more force than (A — 1) par- 
ticles on the long (shallow) side (i.e. is a < 1/N)1 If 
a < 1/N, the object will localize at a position where 
none of the particles is on the steep part of the poten- 
tial. If a > 1/A^, the force exerted on particles in the 
shallow part of the potential will push the localization 
position further to the +x direction. (2) Is the separa- 
tion between particles longer than the short side of the 
sawtooth potential (i.e. is d > aL)7 When d < aL, 
each period of the effective potential U{xcm) will have 
only one minimum. For d > aL, there are A minima 
for each period of U{xcm), each corresponding to one of 
the particles localizing at the minimum of the sawtooth 
potential V{x). 
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FIG. 2: (a) Schematic of ID ratchet potential V{x) with 
a = 0.25 (top) and the effective potential U{xcm) of a dimer 
(d = 0.2L) in this potential. The dimer's localization position 
is x-mi-n = (1 — a)L — d/2. The diffusion distances necessary 
for localization in the adjacent wells in the -l-a; direction and 
in the —x direction are labeled as Aa;+ and Ax- respectively, 
(b) V{x) and U{xcm) for a = 0.25 and a dimer of length 
d = 0.33L. Since now d > aL, there are two possible localiza- 
tion positions. The dimer's average localization position falls 
somewhere in the shaded region (see text). 



Based on these factors, we discuss the behavior of a 
trimer in the following four regimes: (I) a < 1/N =1/3 
and d < aL; (II) a < 1/7V = 1/3 and d > aL; (III) 
a > 1/N = 1/3 and d < aL; and (IV) a > 1/N = 1/3 
and d > aL. 

Regime I (a < 1/3, d < aL): The force on one 
particle on the steep side of the potential is stronger 
than the force on two particles on the shallow side of 
the potential, and therefore the trimer will localize with 
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FIG. 3: Brownian dynamics simulations, (a) Average velocity 
as a function of a for a monomer and for dimers of length d = 
L/5, L/4 and L/3, using kT/Vm^^ = 1/50, L = 5.6 and tan = 
tofi = 20. (b) Critical asymmetry, Oc, as a function of d/L. 
Points below d/L — 1/4 are compared with the prediction of 
Eq. ©, Qc = 1/2 - d/L (solid line). 



one of the particles at the sawtooth potential minimum 
(x = (1 ~ a)L). Because the spacing between parti- 
cles is not large enough for the trimer to straddle ad- 
jacent potential wells, the only center-of-mass localiza- 
tion position is a;inin = (1 ^ Oi)L — d (see Fig. 01a)), 
and therefore Aa;+ = L + d — -Xmin ~ otL + 2d and 
Aa;_ = Xmm — d — (1 — a)L — 2d. The condition 
Ax+ = Ax- yields 
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FIG. 4: Trimers constrained to ID motion (fcT <C Vmax) are 
shown at their localization positions in a ratchet potential 
V{x) for the cases: (a) a < 1/3 and d < aL; (b) a < 1/3 
and d > aL; (c) a > 1/3 and d < aL; and (d) a > 1/3 
and d > aL. Dashed line indicates how the center-of-mass 
position of the trimer corresponds to a local minimum of the 
effective potential U{xcm)- 



Regime II {a < 1/3, d > aL): Since the spacing be- 
tween particles is greater than aL, the trimer now has 
three possible localization positions (Fig. El[b)). The 
asymmetry condition for reversal depends on the rela- 
tive probability for the trimer to localize in each of these 
positions. For a < 1/2, the probability is greatest for 
the trimer to localize in the position furthest in the —x 
direction (a;,nin = (1 ~ ck)L — d). Therefore, it is reason- 
able to expect Ax+ > Ax- for some a < 1/2, yielding a 
critical asymmetry ac < 1/2. 

Regime III (a > 1/3, d < aL): The force on two 
particles on the shallow side of the potential is now 
stronger than the force on one particle on the steep 
side of the potential, leading to the localization position: 
a^min = {l-a)L (Fig.Uc)). Now, Ax+ = L-Xmin = aL 
and Ax- — Xaiin = (1 — Oi)L, which yields ac ~ 1/2, so 
the velocity is always positive in this regime. 

Regime IV (a > 1/3, d > aL): A trimer with d < 
xcM < L + d at the end of toff will localize at x,nin = 
{2—a)L—d, such that the leading bead is at the minimum 
of a V{x) potential well (Fig.gjd)). Now, Ax+ = L + d- 
Xmin = 2d— {l — a)L and Ax_ = x^in — d = {2 — a)L — 2d, 
yielding ac — 3/2 — 2d/L. Since we have d/L > a> 1/3, 
the reversal condition is never met for this regime, and 
the velocity is always positive. 

These predictions can be summarized as follows: The 
trimer has positive velocity for 1/3 < a < 1/2. When 
a < 1/3, ac decreases with increasing d/L. When d/L 
is large enough that ac < 1/3, then the trimer velocity 
reverses direction twice as a is increased from zero to 
1/2: first, from positive to negative velocity at a = ac, 



and second, from negative back to positive velocity at 
a = 1/3. However, if d/L is small enough that ac > 1/3, 
then the trimer does not exhibit a velocity reversal with 
a for a < 1/2, and the qualitative behavior of a monomer 
is recovered. Simulation results shown in Fig. [S] confirm 
these analytical predictions. 

We have shown that the effective potential for a dimer 
in the presence of a 'sawtooth' potential can have up 
to two minima in each period, and it is possible for the 
velocity to reverse once in the range < a < 1/2. Like- 
wise, a trimer can have up to three minima in its effective 
potential and, for the right choice of parameters, the ve- 
locity reverses twice in the range < a < 1/2. In Fig. |^ 
simulation results show that a chain with iV = 4 can un- 
dergo three velocity reversals inO<Q;<l/2, and a chain 
with N — 5 can reverse direction four times in this range. 
In general, if we consider a ID chain of N particles with 
{N — l)d < L, each period of the effective potential can 
have up to N minima, and we expect (N — 1) possible 
velocity reversals. 
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FIG. 5: Brownian dynamics simulations. Average velocity 
as a function of a for trimers of length d — 0.051/, 0.15L, 
0.25L, 0.35L and 0.45L, using fcT/Vmax = 1/50, L = 5 and 
TON = TOFF = 20. The dashed-line interpolation between 
data points is included as a guide to the eye. 



IV. TEMPERATURE DEPENDENT COUPLING 
EFFECTS 



Here, we discuss the temperature dependence of flash- 
ing ratchet velocity for mechanically coupled particles in 
ID, compared with that of individual point particles. As 
we will show, the velocity of a dimer reverses direction 
twice as temperature increases, while a monomer in the 
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FIG. 6; Brownian dynamics simulations, (a) Average velocity 
as a function of a is plotted for a chain of A'^ = 4 particles, us- 
ing d = 0.2L, kT/Vm^y, = 1/50, L = 5 and tqn = roFF ^ 20. 
The dashed-line interpolation between data points is included 
as a guide to the eye. (b) Average velocity as a function of 
a is plotted for a chain of A = 5 particles, using the same 
parameters as in (a). 



same system does not reverse with temperature. First, 
we propose a reason for the ratchet velocity of a dimer to 
reverse direction with temperature, and then we present 
BD results demonstrating this behavior. 

The velocity of an object in a flashing ratchet depends 
on the temperature of the surrounding heat bath in two 
ways: (1) Changing the temperature affects the magni- 
tude of ratchet velocity by changing the diffusion con- 
stant, D = kT/Njb, of the chain. (2) An object in V{x) 
localizes exactly at the minimum of U{xcm) only for the 
limit kT <^ Knax- As temperature increases, the ob- 
ject's Boltzmann probability distribution broadens and 
the mean position shifts away from the effective poten- 
tial minimum, Xmin, as illustrated in Fig.EI For coupled 
particles, a period of U {xcm ) has several regions of dif- 
ferent relative slope, and the symmetry of the positional 
probability distribution shifts with increasing tempera- 
ture. This can lead to a change in the direction of ratchet 
velocity. 

For a dimer with x^-m < L/2, where xqm = des- 
ignates an absolute maximum of U{xcm), we expect 
ratchet velocity to be in the —x direction in the limit 
{kT <^ Vmax), as discussed in section ITlI Al However, as 
T increases, we expect two velocity reversals: (A) From 
negative to positive velocity and (B) from positive back 
to negative velocity, for the following reasons: 

Reversal (A): For a dimer exposed to V{x), U{xcm) 
has a shallower slope immediately to the right (-l-a; di- 
rection) of the absolute minimum than to the left (—x 



direction). Therefore, as T increases from zero, the mean 
localization position increases, and a finite temperature 
can be chosen such that the dimer is more likely to be 
in the region (L — A < xqm < L) during ton than 
in the region (0 < xqm < A), where A is the dimer's 
average diffusion distance during ioff, and is given by 
A = y/2D toff. In this case, ratchet velocity is in the -l-a; 
direction, indicating a reversal from negative to positive 
velocity for increasing temperature. 

Reversal (B): Because a period of U{xcm) has sev- 
eral regions of different slope, continuing to increase the 
temperature can shift the probability distribution such 
that the dimer is more likely to be found in the region 
(0 < xcM < A) than in (L — A < xcm < L), producing 
a second reversal from positive back to negative velocity. 

The physical reason for these reversals is that, at dif- 
ferent temperatures, particles will sample regions of the 
potential that have different slopes. Because the effective 
potential for coupled particles is more complicated than 
the simple sawtooth potential V{x), changing the tem- 
perature can change the symmetry of the center-of-mass 
probability distribution, producing a reversal in ratchet 
velocity. 





U(Xcm) 
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FIG. 7: Upper panels: A ratchet potential, V(x), with 
a = 0.35 is shown (solid line) along with the Boltzmann po- 
sitional probability distribution (dashed line) for a monomer 
exposed to V{x). Lower panels: U{xcm) is shown with the 
corresponding probability distribution for a dimer of length 
d — 0.2L. Probability distributions displayed for two choices 
of temperature: (a) fcT/Vmax = 1/50 and (b) kT /V,-n^^ = 2. 



This prediction of two velocity reversals as a function of 
temperature for a dimer is confirmed by BD simulations, 
shown in Fig. |S1 along with the ratchet velocity of a single 
particle as a function of temperature. As temperature 
increases from T — Q, the dimer undergoes a reversal 
from negative to positive velocity at a point labeled A. 
By further increasing the temperature, another velocity 
reversal is induced, at point B, from positive back to 
negative velocity. By contrast, the monomer always has 
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FIG. 8: Brownian dynamics simulations. Average velocity as 
a function of temperature T for a monomer and for a dimer of 
length d — L/5, using a — 0.35, L = 5.6 and ton = toff ~ 20. 
The dashed-line interpolation between data points is included 
as a guide to the eye. Current reversals are labeled as points 
A and B (see text). 



positive velocity for a < 1/2. 



V. COUPLED PARTICLES IN THREE 
DIMENSIONS 



is straightforward to show that the net linear force on 
a uniform shell with diameter d < L exposed to V{x) 
is zero if a fraction a of the diameter is to the left 
of a potential maximum. The probability distribution 
of a freely rotating dimer meets this requirement if it 
has xcM ~ d/2 — ad or xqm — L + d/2 — ad. Thus, 
Aa;+ ~ L + d/2 — ad — x^i^ aL + d/2 — ad and 
Aa;_ = a;,„i„ - {d/2 - ad) = (1 - a)L - d/2 + ad. The 
symmetry condition, Aa;+ = Aa;_, yields ac = 1/2, so 
there is no velocity reversal for a < 1/2. 




FIG. 9: A freely rotating dimer {d < aL) at fcT/Kiax < 1 
experiences on average no net linear force if xcm = d/2 — ad 
of a;cM — L + d/2 — ad. The dashed-line circles indicate 
the positional probability distribution of a dimer centered on 
these Xcm positions at the end of tog. The average diffusion 
distances necessary to localize in an adjacent well are labeled 
as Ax- and Ax+. 



Here, we examine the behavior of a dimer in a 3D 
system in the limit {kT ^ Vmax), in the same ID external 
ratchet potential, V{x), as in previous sections. When 
a dimer is allowed to rotate freely in 3D, the effective 
potential depends on the dimer's orientation with respect 
to the X-axis. The direction of velocity can be predicted 
based on the following insights: (1) When V{x) is applied 
to a dimer with a random initial distribution, the dimer 
experiences both a linear force and a torque. In response 
to these forces, both particles in the dimer will localize at 
exactly the minimum of V{x), and therefore the center- 
of-mass localization position is Xmin = (1 ~ a)L. (2) If 
tos is long enough that the dimer's orientation at the 
end of toff is uncorrelated with its orientation at the 
beginning, the average particle distribution of a dimer at 
the end of toff is a uniform shell with diameter d. 

Based on these observations, it is possible to calcu- 
late the diffusion distances, Ax+ and Aa:_, necessary 
for the dimer to localize in the adjacent well in the +x 
and ~x direction, respectively, for the majority of possi- 
ble dimer orientations at the beginning of ton- To find 
these distances, we must determine the position xqm for 
which the net force on the dimer, averaged over pos- 
sible orientations, is zero in the presence of V{x). It 



This prediction is confirmed by BD simulations 
fFig. llO|l . Results from Fig. O are included in this plot to 
illustrate how the behavior of the freely-rotating dimer 
matches that of the monomer, in contrast to the velocity 
reversal that takes place for a dimer in ID. We expect 
that the velocity reversal observed in the ID system will 
also vanish for rigid chains of three or more particles in 
a 3D system, because the symmetry arguments used for 
the freely-rotating dimer will also apply to longer chains. 



VI. CONCLUDING REMARKS 

We have shown that, in ID, for kT <C T^nax, a rigid 
chain of evenly spaced particles in a flashing ratchet re- 
verses direction multiple times as a function of chain size 
or ratchet asymmetry. The physical reason is that cou- 
pled particles in a simple sawtooth potential arc effec- 
tively equivalent to a single particle in a more compli- 
cated potential. In this sense, our results are related to 
the finding that a single particle in a ratchet driven by 
dichotomous force fluctuations can undergo multiple re- 
versals when the external potential has multiple wells in 
each period . 
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FIG. 10; Brownian dynamics simulations. Center-ol-mass ve- 
locity as a function of a for a freely-rotating dimer of length 
d = 0.2L, using fcT/Vmax = 1/50, L = 5.6 and ton = toff = 20. 
For comparison, data points from Fig. 3a for a monomer and 
for a dimer confined to one-dimensional motion are also dis- 
played. 



Because a period of the effective potential of a dimer 
in V{x) has multiple regions with different relative slope, 
the center-of-mass probability distribution can reverse 
symmetry as a function of temperature. The shifting of 
the probability distribution during ion changes the likeli- 
hood for the dimer to diffuse to an adjacent well in either 
direction during iofi- For this reason, a dimer can un- 
dergo multiple velocity reversals with temperature. Note 
that these reversals happen for a qualitatively different 
reason than the reversals described for coupled particles 
in the low temperature limit. In that limit, the direction 
of velocity is determined by the symmetry of the location 
of the confining minimum of the effective potential, rel- 
ative to the absolute potential maximum. As a function 
of temperature, on the other hand, reversals occur be- 
cause of the symmetry of the steepness of the confining 
potential. Due to the complexity of effective potentials 
for iV > 1, the two reversal mechanisms are not always 
directly related. 

While the velocity reversals observed in a ID system do 
not occur for freely rotating rods in 3D, this does not rule 
out the possibility for reversal when the object's 3D rota- 
tion is partially confined. The ratchet transport of rigid 
rods confined to a tunnel with radius on the order of the 
rod length could provide a model for coupling effects that 



can be experimentally realized more readily than a true 
ID system. For example, transport of single-stranded 
DNA fragments has been observed using a 2D array of 
asymmetrically spaced, micron-scale electrodes to create 
a flashing ratchet on a silicon chip [s^l- A similar tech- 
nique may be used to experimentally realize the present 
theory. For example, one may expose a polyelcctrolyte 
that is confined in a quasi-lD nanochannel ^5, 36, 3^ 
to a time-dependent, asymmetric electrostatic potential. 
Such a system could provide a method to separate parti- 
cles of different sizes in opposite directions. 

Our investigation of the fiashing ratchet transport of 
mechanically coupled particles may also be relevant to 
the directed intracellular transport of linear molecular 
motors 13, ^3§\ . It has been observed that structurally 
similar molecular motors can move in opposite directions 
m ES El m El El- Binding between a motor head 
and its track is usually by multiple bonds, which one may 
choose to model by several, semi-rigidly interconnected 
"beads" exposed to a binding potential. In dimeric mo- 
tors, an unbound motor head is not free to rotate in 3D 
because of its attachment to the track by another mo- 
tor head. The binding potential is thus quasi-lD, as dis- 
cussed in the present paper. Our observation that objects 
with more internal structure have a larger number of ve- 
locity reversals may explain why, in complex motors, a 
relatively small change in structure can produce reversal 
of walking direction. 

The current study deals with rigid rods which are 
shorter than a spatial period of the ratchet potential. In 
another recent study, we explored the role of fiexibility 
in fiashing ratchet transport by modeling a 3D Brown- 
ian motor based on a polyelcctrolyte or polymer carrier 
with a radius of gyration on the order of several ratchet 
periods, demonstrating that flexibility can increase the 
speed and the stall force of the motor |3- The unique 
collective effects observed in each of these cases suggest 
that it would be interesting to investigate an interme- 
diate regime, such as self-avoiding, flexible chains in 3D 
with total contour length less than a ratchet period. 
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